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' A method of resummation of infinite series of perturbation theory diagrams is applied for studying 

i the properties of random band matrices. The topological classification of Feynman diagrams, which 

was actively used in last years for matrix model regularization of 2d-gravity, turns out to be very 
useful for band matrices. The critical behavior at the edge of spectrum and the asymptotics of energy 
level correlation function are considered. This correlation function together with the hypothesis 
about universality of spectral correlations allows to estimate easily the localization length for eigen- 
vectors. A smoothed two-point correlation function of local density of states p(Ei,i)p(E2,j) c as 
well as the energy level correlation for finite size band matrices are also found. As d-dimensional 
generalization of band matrices lattice Hamiltonians with long-range random hopping are considered 
as well. 

Q ■ PACS numbers: 05.45. +b, 72.15.Rn 
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I. INTRODUCTION 



Random band matrices were introduced many years ago by E.Wigner 0] as a model Hamiltonian for complicated 
quantum systems. In the last few years statistical properties of random band matrices have become again the subject 



of intensive analytical and numerical investigation |2|-|4| due to their application to condensed matter physics and 
I , statistics of spectrum of chaotic systems. 

Up to now all the analytical results for this quasi-ld systems (for review see were obtained by mapping them 
onto a super-symmetric cr-model ||. However, in this paper we would like to develop another method for calculations 
with random banded matrix ensembles. Roughly speaking our method consists of summation of infinite series of 
i— i" perturbation theory diagrams. Diagrammatic methods were used many years ago for investigation of Gaussian 
ensembles of N x N matrices but later this approach was almost forgotten for years. Our aim in this paper will be 
' to show how this "ancient" method may lead rather easily to new results for band matrices. 

Let us consider a Gaussian ensemble of random band matrices. Due to the Wick theorem such ensembles may be 
■ completely defined by the second moment: 

O ! 

O ■ HijHkn = F(i - j)(5j k Si n + QS jn Sik) ■ (1) 

\Q ' The function F(i — j) = F(\i — j\) vanishes very rapidly outside the band (at \i — j\ > b ^> 1). Parameter takes 
0^ ■ values or 1. If = 0, one is dealing with Hermitian matrices of general form (Gaussian unitary ensemble - GUE[3), 
while = 1 corresponds to the real symmetric matrices (Gaussian orthogonal ensemble - GOE). It is convenient to 
define the width of the band b and the typical strength of the interaction V through moments of the function F: 



£ ■ .Fo L,F{n) 

g : V^F Q = Y^F(n 



For practical computations we will sometimes use F of the form 



As will be shown below the results essentially do not depend on the details of the shape of function F(i — j). We need 
only F(i — j) to be sufficiently smooth so that after averaging (0) all discrete sums may be replaced by integrals up 



1 Although, the notations GUE and GOE for our ensembles are mainly the traditional since unitary (orthogonal) invariance is 
broken explicitly for banded matrices. 
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to negligible corrections ~ exp(— b). In doing so we still are able to consider the corrections of any finite order in 1/6. 
Moreover, we argue below that even smoothness of F does not seem to be necessary for most interesting applications 
(see ( |l5| , p^| ) and discussion below, it will be shown also why (||) is the most natural choice of F). 

By simple d-dimensional extension of band matrices one obtains the Hamiltonian for particle hopping on a d- 
dimcnsional lattice with random nonlocal interaction. This lattice model may be described by the same formulas 
(0H) with simple replacement of all integer indices by the integer d-vectors and trivial redefinition of F 

^-^ = ^ eX p("w) ' (4) 

Contrary to band matrices, properties of this model seem to be completely unknown. The b 3> 1 in (||) is effectively 
the number of "neighbors", connected to each lattice site. The band matrices ([!]) may be associated now with the 
random Hamiltonian for one dimensional lattice. 

The "physical quantities" which we would like to consider are connected with the Green's function and the local 
density of states 

GijiE) = (^Th) , p{E,i) - hmG^(E - iO) . (5) 

There is no summation over i in Gi^ in the last formula. More specifically we would like to consider the averaged 
density of states p{E) and the correlation of densities for different but very close energies (and even for different 
positions i and j). 

It is also constructive to compare our results with those for the ensembles of usual N x N random matrices, which 
are defined by the second moment 



_ V 2 V 2 

F(i- j) = — ; HijH mn = — (S jm S in + QS jn S im ) . (6) 



This Hamiltonian may be considered as the d = reduction of the lattice model (Q) . Historically three main 
approaches were applied for studying the statistics of the full N x N matrices. The description of this approaches may 
be found e.g. in papers The first one is summation of infinite power series in V/E Q. Another two methods 

are the replica trick J§| and the super-symmetry method || . 

The success of both replicas and super-symmetry is essentially based on the use of Hubbard-Stratonovich trans- 
formation. For N x N matrices this transformation reduces the problem to an almost trivial calculation of a few 
dimensional integral. On the other hand for band matrices even after Hubbard-Stratonovich one still is faced with a 
cr-model on one-dimensional lattice. Thus it seems quite probable that neither replicas, nor super-symmetry will lead 
to considerable progress in the lattice model fl). 

The exact solution of 2d quantum gravity [ flO| stimulated the explosion of interest in matrix models. In this appli- 
cation of random matrices the discretized random surfaces appear as Feynman graphs in the perturbative expansion 
of the matrix integral. However, technically the famous double scaling solution of 2<i-gravity has nothing to do with 
summation of graphs. For ensembles invariant under orthogonal transformations it is useful to work with N eigen- 
values instead of all N 2 matrix elements. Unfortunately, for band matrices, or lattice hopping Hamiltonian we could 
not find such a simple solution, which is not based on the diagrammatic expansion. Nevertheless for models ([j], f|) 
the topological classification of diagrams, which arose in 2<i-gravity (and originally in QCD simplifies drastically 

the summation of the series. 

This experience of dealing with the matrix integrals for 2<i-gravity was later used for solving the problems typical 
for quantum chaos. In [Q a method of calculation of correlators of Green functions for ensembles of large N x N 
matrices was developed. An approach based on summation of perturbation theory series for various ensembles of 
random matrices was also used in the series of papers of Brezin and Zee (see e.g. [|l3| and the later papers of the same 
authors), though their diagrammatic technique differs from used in the present paper. 

The statistics of band matrices with parameters of the band slowly varying along the diagonal was considered in 
fL4f . The behavior of the edge of spectrum for this extensions of the model (Q) showed some surprising similarities with 
the edge properties of matrix ensembles considered for the 2<i-gravity. The topological classification of the diagrams 



which we will explore below was also briefly discussed in 14 



The organization of this article is as follows. The general description of the diagrammatic technique is given in 
Section 2. By comparison with N x N case a diagrammatic proof of semicircle density of states is found. We also 
develop partial summation of infinite subseries of topologically trivial tree-type diagrams. In Section 3 the ideology of 
double scaling limit is used to study the edge of spectrum for random band matrices. The edge of spectrum for lattice 
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model (H|) is considered in Section 4. Surprisingly the critical behavior at the edge for lattices with random hopping 
coincides with the critical behavior of the string-theory inspired model considered in 16 1. In Section 5 the two-point 
correlation function p(Ei)p(E 2 ) c is calculated. More precisely we found the so called smoothed correlation function 
in the large b limit and the first ^1/6 correction to it. Moreover together with the hypothesis about the universality 
of spectrum fluctuations Jl^] this correlation function allows to find the correct estimate of the localization length for 
the eigenfunctions of band matrix. This universality holds also for the ~ 1/6 correction to the correlation function, 
even though the ~ 1/6 correction itself turns out to be a subject of strong cancellations. Finally in the Section 6 
some quantities which have no analog for usual N x N matrices arc considered. These arc the correlation function 
of local density states p{E\, i)p(E2, j) c and the usual density-density correlation function for the banded matrices of 
finite size N. 



II. DIAGRAMMATIC TECHNIQUE 

It follows immediately from ([!]), that only diagonal terms survive in the averaged Green function 

G~ j =G{E)8 itj . (7) 

Let us expand G(E) in a formal series 

G = I^3T = _L YfrfB-Y . (8) 
N E - H NE ^ \E y ' 

n=0 v 7 

Now all we need is using the Wick theorem and the second moment (|l|) to calculate all the average values of the 
trace of product of H - matrices. In the standard Feynman diagram technique each H n corresponds to n-legs vertex 
and the averaging in (||) reduces to counting of the number of possible ways of contracting these legs to each other. 
However, like in 2d-gravity it is more convenient to draw the dual Feynman graphs. For dual diagrams each ify 
corresponds to a segment with numbers i and j at the ends, while H n corresponds to n-vertex polygon with matrix 
indices i\, 12, ■ ■ ■ i n assigned to the vertices (see fig. la). It is also useful to draw the arrow on each segment showing 
the direction from first to second index. Within this language the Wick contractions in (||) correspond to the glueing 
of pairs of segments. Our aim now is to calculate the number of ways in which the edges of the polygon may be glued 
into some closed surface. 

For Hermitian matrices (O = in ([l])) the segments should be glued in the opposite direction thus forming the 
oriented surface. For symmetric matrices (9 = 1) the nonoriented surfaces are also allowed (e.g. the Mobius band). 





a b 

FIG. 1. One plaquette a. and sphere b. glued from such plaquette 



The example of simplest surface of spherical topology is shown in fig. l.b. It is easy to verify that just the spherical 
surfaces dominate in G at large b. Only in this case the sum over n matrix indices for H n gives the factor Nb° and 
thus (see (|),(|)) 

G spherical rv ' & ■ (9) 
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Moreover, as may be seen from fig. lb, the summation over each index in spherical diagrams is completely independent 
and results in a fixed factor V 2 for any choice of the function F(i — j) (||) and for each of ensembles ([[]) , (0) , ^) . In par- 
ticular this means (analogous observation was also done in ]l3| ) that the Green function in the leading approximation 
coincides for all ensembles (0), 



g = —(e-Ve 2 -w 2 ) (10) 

(a very clear proof of this formula for full N x N matrices see in 

§§)• „ 

Before considering the ~ 1/6 corrections let us carefully examine equation (|l(]). As we have said above, Go may 
be thought as the exact sum of the part of the series (||) corresponding to diagrams of spherical topology. It is seen 



from (10) that this series is convergent only outside the circle on the complex E plane with a radius \E\ = 2V. The 
values of Go inside the circle (\E\ < 2V) may be found via analytic continuation. This feature of the series has two 
important physical consequences. Firstly, if one is approaching the singular points E = ±2V starting from large E 
the more and more complicated diagrams became important thus approaching some kind of continuum limit. It may 



be not so obvious from the equation (10) because for Go the series converges even at E = ±2V but as we will see 
below the ~ 1/6 corrections to G are more singular and the summation of (V/E) n do is saturated by the terms with 
very large n. 

Even more troublesome turns out to be the calculation of correlation functions. In this case one has to consider the 
Green functions close to the borders of the cut (E = E o ±i0 , \E\ < 2V) and very far from the domain of convergency 
of the series^). The price for such unreliable procedure will be the severe cancellations of different contributions to 
the ~ 1/6 corrections to correlation function (see Section 5) 



As we have taken into account exactly in (10) all spherical contributions, the ~ 1/6 corrections naturally turn out 
to be determined by the diagrams of the more complicated topology. As we have seen before @,@ the sum over 
spherical graphs for full N x N matrices and for band matrices coincides up to trivial replacement N — > 6 because 
the summation over each matrix index in the tree diagram of the kind of fig. lb is independent and results exactly in 
the trivial factor V 2 . This direct correspondence does not holds for the diagrams of the more complicated topology. 
However, each sum over matrix index contains effectively ~ 6 items and their magnitude corresponds to that for 
N x N matrices up to substitution 6 — > N . Thus all the business with classification of the diagrams in the powers of 
the parameter 1/6 (1/N) , which was so productive for the full matrices, still holds for the band matrices as well. In 
particular for the Hermitian band matrices (0 = in (hj)) one may use the well known Eulers theorem to show that 
the corrections to ( fl0| ) may be only of the kind ~ (1/6^" with n being the number of handles. 

Up to now we have associated each Feynman diagram with some surface. However, as it may be seen e.g. from the 
fig. 1 because in our problem we have in fact one large plaquette (or two for the correlation functions) it is natural 
to consider only the border of this plaquette which should be glued to some kind of branched polymer. It is seen 
from the fig. lb that the spherical surfaces are associated with the tree-type polymers. On the other hand the ~ 1/6 
corrections to G will be associated with selfintersections (the closed loops) of the polymer. It seems very attractive 
to divide the calculation of the ~ 1/6 corrections into two stages. The first stage will consists in summation over the 
trees and at the second stage one will take into account only the dressed selfintersected diagrams. To this end it turns 
out to be useful to consider instead of the Green function the logarithm 



i* h(i-f) , ow-i + §- (ID 

The simple combinatorial calculation allows to replace the perturbative series for L{E) by the sum over skeleton 
graphs, as it is demonstrated in the fig. 2 



n \E J N ^ p 

n ~ no contractions 

of neighbours 



Here the contractions between the nearest neighbours in the last sum are forbidden because they have been taken 
into account exactly in Go-s. More precisely for the real symmetric matrices (0 = 1 in (Q)) the contractions between 
neighbours still are allowed but only due to the second term (~ O) in the r.h.s. of ([!]). 
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FIG. 2. Reduction of ln(l — E/H) to skeleton diagram. 



III. THE EDGE OF SPECTRUM 



Now at last we are able to consider the edge behavior of the ~ 1/6 corrections to the Green function Go (pOj). To 
this end in particular we have to take into account the long chains of glued dressed links of the kind of ([l2]) (or fig. 2). 
Consider the simplest two-link chain which is shown in the fig. 3: 



Here in order to calculate ^ n we have used the specific form of the function F(i — j) (||). The various methods may 
be used in order to prove (|l3|). For example one may use the mathematical induction method. While deriving the 
equation ( |l3| ) we have replaced all the summations over the intermediate indices z x , z 2 , ... ,i n -i by the integrations. 
The accuracy of such procedure S 1 ^ ~ e~ b for any smooth function F(t — j) still allows one to consider the ~ (l/&) n 
corrections for any n < b . 



^n{in — io) — ^ {Hi a i 1 Hi 1 i )(H il i 2 Hi 2 i 1 ) . . . {H in _ lin H in i n _ 1 ) 



(13) 



Jl22...«n-1 





FIG. 3. Double-line chain calculated in (13). 



The ^ n evidently satisfies the sum rule 




(14) 
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which also may be used in order to find the normalization of (|13[). 

The sum rule ( |l4|) holds exactly for any choice of the function F(\i—j\) while the equation jis] ) is model dependent. 
However, as it will be shown now, for large n the equation (|3|) is also universal. Consider the recursion formula for 
large n (and arbitrary F) 

y n+1 ( i ) = Y,F(\ i -3\)y n ( J ) = V 2 y n (z) + —^ n (i)+ ... . (15) 

3 

This equation is the discrete (in time) analog of the heat conductivity equation. Together with the initial condition 

*o(i)~*(0 , (16) 

and the sum rule jl4| ) equation ( |l5| ) allows to reproduce the formula (|l^) for the chain. In fact this is the reason for 
considering the function F(\i — j\) of the form (|3|) as the most universal one. Moreover, even if one starts with some 
irregular function F (which may seem to be crucial because only for smooth F(x) the summation may be replaced 
by integration with the accuracy ~ e~ b ), taking into account of a long chains (n 3> 1) effectively smooth it out. 

Up to now we have considered only the band matrices. However the formula ( |l3| ) may be easily generalized for the 
random-bond lattice case (Q) 

V 2n f 1 r2 



*^ )= R<^ exp r^ > ■ (i7) 



FIG. 4. The Mobius band corresponding to the ~ 1/6 correction to Green function for real symmetric matrices. 



At this point we have completed all the preliminary formalities and are able to calculate the ~ 1/6 correction to the 
Green function. Let us consider the ensemble of real symmetric matrices(GOE) for which the correction of the first 
order ~ 1/6 exists. The corresponding skeleton Feynman graph is shown in the fig. 4. It is seen from the figure that 
the collinear links should be glued in the (1/6) 1 order and thus this diagram is forbidden for the Hermitian(GUE) 
matrices. In terms of surfaces fig. 4 corresponds to the Mobius band. Combining together (11), ( |l2| ) and (13) one 
finds the correction 



1 c f H 



N ^ 2p 

P =i 



skeleton — 



1 °° 1 



b^2p 

P =i 



(18) 



or for the Green function 



G x = SG = 



n i oo 
Lr \ - 



(G V) 



2p 



b(2n) d / 2 Go ^ P 



d/2 



(19) 



We give the result at once for arbitrary dimensionality < d < oo. For band matrices one has to choose d= 1. 

The only important factor which is responsible for the difference between usual N x N matrices (d = 0), band 
matrices and, random hopping Hamiltonians in ( |l9| ) is the \/p d / 2 in the sum. Technically this factor comes from the 
double chain of matrix elements (O,^). One may consider the length of the Mobius band (fig. 4) p as the discrete 
time in the diffusion equation ([151). Then the l/p d ^ 2 will simply correspond to the time dependance of the return 
probability for classical diffusive particle at the time p for different dimensions. 



G 



We are mostly interested in the singularities of G at GqV — > 1 (or E — > ±2V, see (|ld|)). It is seen from (18) that 
Li has a finite limit at GqV = 1 for any d > 0. On the other hand the G\ (|l^) is convergent and the diagram of fig. 4 
approaches the continuum limit, at least for d = 1 and d = 2 (for d > 2 the G\ is also singular, but mostly due to G' Q 
without any continuum limit). 

For the random band matrix case (d = 1) it is easy to found from (K%, (fL0|) that close to singular point 



G(£7 



2vo = ±-A 



£-2v 



V 



1 1 

6 17 



V 



E-2V 



3/4 



1 const 



b 2 V \E-2V 



(20) 



Here the last term is the order of magnitude estimate of the ~ 1/b 2 contribution. The Feynman diagrams corresponding 
to this contribution are shown in the fig. 5 (the explicit calculation of one of them will be presented below). 






a 



FIG. 5. The three possible second order ~ 1/b 2 diagrams for ensemble of real symmetric (GOE) matrices, a. - dumb-bell 
shaped diagram, c. - torus 

The simple counting of the power of convergency for the higher order diagrams shows that close to singularity the 
Green function should be described by some scaling function 



v vb 2 / 5 V V v 

So, one may conclude that the singularity at the edge E — ±2V of the perturbative Green function ( JlCj ) should be 
smoothed out at the distances AE ~ V6~ 4 / 5 from the singular points. For example it is easy to estimate the number 
of energy levels falling into this region. The same estimate evidently holds for the number of the levels outside the 
circle \E\ > 2V 

AN - NAE 3 / 2 - . (22) 

This estimate is of particular interest because the density of eigenvalues outside of the circle is purely nonperturbative 
and could not be found in any finite order over 1/b. It is to be noted that for usual N x N matrices (both GOE and 
GUE) AE - iV" 2 / 3 and AiV - 1. 

As we have told in the introduction, random band matrix may be considered as a Hamiltonian for particle on Id 
lattice with random hopping. Eigenvectors for such Hamiltonian should naturally have some finite localization length. 
This localization length for the band matrix ensembles (jlj) was found in the papers of Fyodorov and Mirlin within 
the super-symmetry method (the calculation of li oc in our approach will be given in Section 5) 

hoc ~b 2 (l- ^) , \E\ < 2V . (23) 



This result in fact was found only in the leading order in 1 /b and should be changed at the edge of spectrum. One 
may combine our result (pi]) with ( p3| ) in order to estimate the localization length for nonperturbative states living 
outside of the main band (circle) 
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C ~ b 2 AE ~ 6 6 / 5 

IOC 



(24) 



In particular l^ c AN ~ AT which may mean that nonperturbative states spatially do not overlap. 

In the recent paper fL5| the distribution of Lyapunov exponents for random band matrices was studied numerically 
at the very edge of perturbative part of spectrum \E\ — 2V. The authors of ref. have considered the band matrices 
with step form of the band F(\i — j\ > b) = for which there exist exactly 2b + 1 Lyapunov exponents. Their result 
reads 

A^ft-V 3 * 1 ^ . (25) 

Here i = 1, 2, 3... , and Ai is the smallest decrement of the Lyapunov exponent. The corresponding solution grows like 
ip$ ~ exp(nAi). Naively one may expect that the smallest Lyapunov exponent directly gives the localization length 

Ar x «Jioc • (26) 

However, this simple guess puts in conflict the two results (|24|) and (|2jj). On the other hand, the numerical accuracy 
of || allows one to be sure in (pa) for sufficiently large i, namely 1 <C i <C b. The accurate result for smallest 
Lyapunov exponent still may differ from (p5|). 

Moreover, the direct correspondence between the localization length and first Lyapunov exponent naturally take 
place in the central part of spectrum \E\ < 2V but at the edge \E\ = 2V formula ( p6| ) is not necessary correct. It 
seems rather probable that the eigenstates outside the circle \E\ > 2V appears due to some very rare fluctuations of 
our random banded Hamiltonian Hij. Furthermore, we have no idea how even to estimate the density of such rare 
fluctuations. If the number of this fluctuations is small compared to N/b & / b the eigenstates with \E\ > 2V will not be 
distributed homogeneously, but will be concentrated in rather rare and dense bunches. On the other hand, the first 
Lyapunov exponent Ai evidently comes from the whole range of variation of vector index. Thus the localization length 
at \E\ = 2V should not be necessary of the same order of magnitude with the first Lyapunov exponent. Unfortunately, 
if the states with \E\ > 2V do are concentrated together one should most naturally expect that A^f 1 » l^ c . Thus, 
the discrepancy between (pll) and (E9) turns out to be even severe in this case. 




i kj j 

FIG. 6. Gluening of dressed hexagon for ~ 1/b 2 correction for Hermitian matrices (torus). 

In fact in (^), ( p~l|) the estimate for higher order corrections was declared without any proof. However, both for 
the more rigorous proof of ( pl| ) and for future calculations it is useful to calculate explicitly at least one nontrivial 
(beyond the one loop) diagram. Therefore we would like to found now the leading ~ 1/b 2 correction for the GUE 
band matrices. In this case the only diagram of fig. 5c survive which corresponds to torus in the surface language. 
The fig. 6 shows how one have to glue the dressed hexagon (see (|l2|)) in order to built this diagram. Indices ki, &2, &3 
on the figure are the lengths of the double-link chains (|l3|), while i and j stands for the matrix indices corresponding 
to the ends of these segments. It is to be noted that both indices i and j appear three times in the diagram, which in 
particular leads to loss of two powers of large b for this diagram compared to the tree one (fig. lb). More precisely, as 
it may be seen also from the diagram fig. I b not the multiple indices but the closed loops are the direct source of the 
1/b . One of the simplest ways to find the order of the diagram in the (l/b) n classification is to look for the number 
of links which are enough to break in order to get the tree-kind simply connected diagram. Now one may combine 
( |l2l ) and (|l^) in order to found the contribution of the fig. 5c 

L " = -Jt E m A 1 1 ^ + fc2 + fca)^(Go) 2fcl+2fc2+2fc3 (27) 
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Here (fci + fc 2 + fcj) -1 comes from 1/p in ( |l2| ) and the combinatorial factor (fci + fc 2 + fc3) +1 takes into account the 
number of positions available, say, for the left end i of the upper segment k\ on the fig. 6 and the right end j of the 
same down segment (these points are shown by the two small circles on the figure). Because these two circles i and 
j are equivalent the combinatorial factor is only ^2 ki , not the 2 ^2 ki . Finally, after one has taken into account by 
this ^2 ki the freedom in definition of the starting point on the circle the sets (fci, fc2, ^3), (fc2, fe, ki) and (fca, fci, /C2) 
became indistinguishable which is taken into account by the factor 1/3 in (p7j). The simple substitution of ( |l3|) into 
@ leads to 

L? = (28) 



J_ v (G y) 2 ^+ 2 ^+ 2fc 3 1 f J_/J_ J_ J_\ 

1 ^ (G y) 2 ^+ 2 ^+ 2fc 3 

Thus the Green function close to the edge of spectrum for the GUE band matrices reads 

2 



1 1 E-2V 1 I f V 



V V V V b 2 V \E — 2V 

1 1 1 iM^)) , (29) 



y F6 2 / 5 V V ^ 

where 

/= f 5{\-x-y-z) dxdydz 

Jx,y : z>o 487r V X V + xz + yz 

In order to be sure that we have done nothing wrong with the combinatorics one may easily repeats the calculation 
( p7| - |29| ) for usual N x N matrices. Some new (unknown) scaling function $ appeared in ([?9j). We have used the same 
argument for $ as for $ ( |2l| ) but the asymptotic series for $(x) is in powers of x~ 5 / 4 , while for in powers of 



IV. THE EDGE OF SPECTRUM FOR LATTICES 

Even more puzzling turns out to be the edge of spectrum behavior for the lattice ensembles (|]). On the one hand, 
the naive power counting for diagrams like, e.g., the diagram of fig. 5c which we have calculated in ( |29| ) gives 

G = I + I^*( 4 Afc2T)) . (31) 



V Vf,th \ \ V 

On the other hand, it is easy to find the first order (~ 1/6) correction to the Green function at the edge from ( |l9| ) 



Gx = —^—\ — - — lnf — - — I , d = 2 , (32) 
1 8nbV\ E-2V \E-2V' 



I V 1 

Gl ^2W\jE~2v]^j2^ ' d>2 • 

It is seen immediately that at least for d > 2 equations ([p|) and ( pl| ) disagree. In fact the solution to this paradox 
was found a few years ago in a paper of Ambj0rn et al. [p. 6| there the toy model for string in d-dimensions has been 
considered. The authors of ref. |l6| have restricted the class of triangulations for string embedded in d-dimensions 
to those having minimal cross-section. While doing so they have obtained effectively the theory of d-dimcnsional 
branched polymer. As we have mentioned above, the summation over dual Feynman diagrams for our band matrices- 
lattices also reduces to summation over some branched polymers. In the lattice case indices assigned to the ends of 



9 



each link expands over d-dimensional (although discretized) Euclidean space while the factor F(i — j) ([|) regulates 
the spatial size of the link just like in the model of lira] . The branched polymers were also considered many times 



within the random vector-matrix models approach (see e.g. 17|). However, only the critical exponents for our model 
of branched polymers (or it is better to say "branched tapes" as it is seen from the figs. 4,5,7) should coincide with 
those for another models. The scaling function itself may be different. 






FIG. 7. Tadpoles or new trees for lattice models. 

In order to solve the contradiction between ( j3l[) and (|32|) it is enough to observe that for d > 2 some of the diagrams 
for the branched polymer are much more singular than others with the same topology. These are the so called tadpole 
diagrams shown in fig. 7. Moreover, each of the diagrams of fig. 7 behaves effectively like some tree diagram (see fig. 
lb) up to some trivial factor associated with the ends of the tree. Therefore, while the true complicated diagrams 
(say of fig. 5b, c) became less and less singular in higher dimensions in accordance with (|3l| ) the tadpole diagrams for 
d > 2 have all the same singularity. On the other hand, one may try to sum up exactly this subseries of rather simple 
diagrams, as it was done in Era. 




FIG. 8. The reduction of tadpole diagram to corresponding correlation function. 

The procedure of calculation of the tadpole contribution is illustrated by the fig. 8. Let us replace each n-tadpole 
diagram by some correlation function of our Green function and (trH) n calculated in the spherical approximation. 
Explicit formula corresponding to such procedure is as follows 
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G 



tadpole 



tr- 



1 (ntrH) n 



E- H 



(33) 



where subscript c stands for connected diagrams. The coefficient k may be found for example from (j32|). The factor 
1/nl in d33| ) accounts for the permutations of various trH. Due to this 1/nl the sum over n results in the trivial 
exponentiation of KtrH. In order to find G 'tadpole it is useful to write explicitly the integration over H 



Gtadpoie = z DH ]y tr E _ h ex P { ~ X! H H H ji M {i - j) + KtrH} 



(34) 



Here M{i — j) ~ F(i — j)^ 1 and Z is the same integral without j^tr E ^_ H but with KtrH included in the argument 
of the exponent. By such choice of Z one gets rid of disconnected diagrams. 

At least in the leading (spherical) approximation formulas ([33]) and (Q) give the same Green function at the edge. 
On the other hand, the Green function ( |34| ) reduces to the zero order one (|l(]) by trivial substitution 



After comparison with (5^) one finds 



Hi 



Hij + const Sij 



E — > E + const 



(35) 



G 




(36) 



d > 2 



For d > 2 this result is almost trivial. Taking into account of the most singular series of corrections result in a simple 
shift of the edge of the cut. For d — 2 taking into account of the tadpoles results not only in the shift of the edge 
by AE ~ ln(6)/6 but also in some nontrivial change of Green function. The singularity in ( |36"| ) for d = 2 will be 
smoothed out at \E — 2V — S ln(6)| ~ ¥■ by some unknown scaling function. 

For GUE matrices (lattices) the simple tadpoles consisting of small Mobius band are forbidden. However one may 
arrange slightly more complicated ~ l/b 2 tadpole shown on the fig. 9 which is allowed also for oriented surfaces. As 
a result for the GUE lattices with random bonds the critical dimension is d = 4 instead of d = 2. 




FIG. 9. The ~ l/b 2 tadpole for Hermitian(GUE) matrices. 



V. CORRELATION FUNCTION 



Physically may be the most interesting quantity which would be calculated with the random matrix ensembles is 
the two point correlation function of density of states at very close energies. It is generally believed that just such 
local quantities most adequately reproduce the measurable features of complicated quantum systems. 

Unfortunately the perturbative procedure, which we are only able to perform, has a serious drawback in the case of 
correlation functions. As it was said in the Section 2, the perturbative series in 1/E are convergent only outside the 
circle \E\ > 2V while the series in l/b turns out to be the asymptotic series. Besides the perturbation theory there 
may exist some nonperturbative contributions say of the form ~ exp(— b(E — 2U) 7 ) with some j ~ 1 . However, after 
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analytic continuation to the border of the cut, which goes from E = —2V to E = 2V , these corrections may be (and 
do are at least for the usual N x N matrix ensembles) converted into some oscillating functions like ~ s'm(b 2 AE) . 
Certainly in our perturbative result all these oscillating terms (if they are) will be smoothed out. For example for the 
usual CUE ensemble one gets instead of the exact result K(AE) = sin 2 (NAE)/AE 2 only K = 0.5/ AE 2 . Due to 
this principal limitation the quantities which we would calculate are called the smoothed correlators. 

The intriguing feature of our diagrammatic series which was considered in the previous sections is that the dressed 
Feynman diagrams at the edge of spectrum approach some kind of continuum limit. Technically it happens because 
the value of (VGq(E)) 2 © which was our actual expansion parameter became equal to (VGo(2V)) 2 = 1 at the 
border. Now, because we are going to work with the energies inside the band |ia 2I < 2V^ there seems to be no room 
for the continuum limit. However, let E\ approach the upper border and E2 approach the lower border of the cut . 
Then it follows immediately from ( [l0| ) that 

VG Q {E X ) = e ^- A / 2 , (37) 
VG (E 2 ) = e -*-V2 ( 



where small A is given by 



TP TP 

A = ■ 1 2 , ReX > . (38) 
i\JW 2 - E 2 



Thus again at least those subseries of diagrams which will have as the expansion parameter the combination 
V 2 Go (Ei) Go (E2) = e~ A may approach the continuum limit (another less direct example of the "effective" continuum 
limit will be considered in the following subsection). 

In order to find the correlation function of two Green functions it is natural again to consider the logarithms 



7 Yin I 1 - — ] Tr In [ 1 " 



Ei J \ Ei 



y^^Tr{{Go{Ei)H)A -Tr{ (G (E 2 )H)p) 

„„/ P I J skelet P I J 



p. v 



(39) 

skelet 



N ^ p 1 



ES^^Go(£i)Go(B 2 )) P 



J b(2n) d / 2 ^ p 2 p d / 2 



p 



Here we have expanded each logarithm in the sum of closed skeleton chains like in eq. (|l2|). Also here and below 
subscript c means "connected" . A factor 2 in front of the last sum accounts for two allowed directions of glued chains 
(Cooperon and diffuson in solid state physics). Finally, one combinatorial factor p in the last sum accounts for the 
number of different ways to contract two skeleton rings of the length p. We see that just the e _A (^8|) turns out to be 
the expansion parameter in (p9t) 



By simple differentiation of (p9j) one finds the correlation function at very close energies E\ — E2 — > 

1 1 9N °° n l - d / 2 

Tr - Tr - = — e~ Ap - (40) 

Ei -H E 2 -H c b(W 2 -E 2 )^ i {2Tr) d / 2 { ' 



Ex- Ei 



2-d/2 



d<4 



This equation reproduces the known result of Al'tshuler and Shklovskii for correlation of energy levels in small 
disordered metallic samples. 

The factor l/p d l 2 in (|3^,[4(t) appears after glueing of two H p into closed two- link chain ([l3],[l7|) . As it was said also 
after eq. ( |l9] ) this factor 1 jprl 2 works effectively as the probability for diffusive particle to return to origin after time 
p. In principle one may go even further in this analogy with classical diffusion. By making the Fourier transform of 
the function ^ n (i — j) ( |l3] , |l7| ) the equation ([lO]) may be written in the form of the integral of squared Green function 
of the diffusion equation (again in agreement with corresponding formulas from jl8||). 

Now let us consider in more details the correlation function for band matrices (d — 1). We would like to show 



how the information contained in the smoothed correlation function (40) combined with the simple hypothesis of 
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universality of spectral correlations allows to find the correct estimate of the localization length for eigen-vectors of 
random band matrices. Consider instead of Green function G the density of eigen-values p(E) , which is simply the 
imaginary part of G (||). It is generally recognized, that the fluctuations of p(E) for all ensembles of full random 
N x N matrices are universal. This means that, being properly normalized, the density of states - density of states 
correlation function has the form: 



P {E l)p{ E 2)c=K r^ {Ei _ _ (4i) 



p(£i)p(£ 2 ) \ dE 

Here the function K(x) ~ 1 at i ~ 1 and dn/dE is the averaged density of eigenstates. Of course, both E\ and E2 
in ( ffl| ) are real (have reached the border of the cut). The detailed form of K(x) is specific for the ensemble under 
consideration (e.g. GUE or GOE), but for given ensemble K is the universal function of energy interval AE = E\ — E2 
measured in units of mean inter-level spacing dE/dN . 

For the random band matrices one may expect the same universal behavior of the correlation of density of state 
fluctuations as (|dj) only if all eigen-vectors are delocalized. On the other hand it is easy to write down the natural 
extension of the formula ( ffl|) for systems with finite localization length I : 



p(E 1 )p(E 2 ) c _ I K (ldn (El _ E2)] . (42) 



p(E 1 )p(E 2 ) N \N dE 

By writing this formula we suppose that the fluctuations of density of energy levels for random band matrices still 
are universal if the energy difference AE = E\ — E 2 is measured in the units of effective mean inter-level spacing 

N dE 

AE eff ~ T - . (43) 



Roughly speaking AE e ff is the mean inter-level spacing for a band matrix of a finite size N ~ I (of course in (41 42) 
we suppose that at least I b). Moreover, if the modified universality do takes place and I ^> b the formulas (42 4^) 
will also account for the energy dependence of the localization length. Of course the two universal functions K(x) in 
(fll]) and (|4|) are completely different. 

Our line of reasoning in fact follows the consideration of (see also |2Cj) for density-density correlator in disordered 
metal. Now we would like simply to convert the arguments of authors of ref. |19| in order to estimate the localization 
length. The averaged density of states for large b may be easily found from (TOT) 

dN 1 _ — N 



p(E) = — = -Im trG(E - iO) = V 'W 2 - E* . (44) 



Now from (|40|) one finds the asymptotics 



p(E 1 )p(E 2 ) c _ -1 V 4 I 



p(E 1 ) p(E 2 ) Nb (41/2 - £ 2 )5/4 \E X - £ 2 |3/ 2 
Together with HWj and (E2L) this correlation function allows to find: 



in accordance with the result of Fyodorov and Mirlin ||] . We have defined here the localization length (or effective 
localization length) by choosing the overall normalization constant in ( f46| ) to be equal to one. Formally, universality- 
based arguments allows to find I only up to some normalization constant ~ 1 , which on the other hand, depends on 
the explicit definition one uses for the localization length. For example, I defined via the inverse participation ratio 
Q or by the first Lyapunov exponent may differ by some trivial factor. Anyway, our estimate 1{E) seems to be much 
less complicated than those of the supersymmetric approach of Q . 



A. Corrections to the correlator 



In this subsection we would like to consider the ~ 1/6 corrections to the correlation function. As it was explained 
above, we are able to calculate only the smoothed correlation functions. More concretely we are able to consider the 
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Green functions not too close to the border of the cut \ImE\ AE e ff ~ V/b 2 ( ^3| , ^6| ). Nevertheless, the uncertainty 
in the correlation function due to smoothening of fast oscillations decrease exponentially like ~ exp(— ImE/ AE e ff). 
Therefore even for smoothed correlation function one is able to consider the corrections of any finite order in 1/6. 

The skeleton diagrams for ~ 1/6 correction to the correlation function (|40|j45| ) are shown in the fig. 10. At a 
first stage one may easily estimate the power of singularity for each diagram at small AE. This singularities are 
associated with the number of summations over the length of tape glued from the segments originated from two 
different logarithms (|9|) ln(l - H/Ex) and m(l - H/E 2 ). Thus compared to the zeroth order result (po| , ^5| ) the 
diagrams of fig. 10a. are of the relative amplitude ~ b~ 1 AE~ 2+d / 2 while the diagrams 10b. are of the relative 
amplitude ~ b~ 1 AE~ 1 . 







a b 

FIG. 10. The ~ 1/6 corrections to correlation function, a. corrections singular like 1/AE 3 for d = 1, b. corrections singular 
like 1/AE 5/2 for d= 1. 



On the other hand, the hypothesis of universality fl4lq ) which allows us to find so successfully the localization length 
(|t|) leads to a strong restriction on a possible form of corrections to the correlation function. As we have considered 
in the Section 2, the spherical result for the correlation function ([l5|) is the exact result in the large 6 limit. Therefore 
the corrections to ( [45|) should be of the relative order ~ 1/6. In its turn due to the universalit y (j42|) 6 may appear 
in the result only in the combination (Ei - E 2 )/AE eSf ~ {Ei - E 2 )l ~ {E x - E 2 )b 2 (see (||g^6|)). Thus the only 
form of the correlation function consistent with the universality condition (|42|) is 



Tr 



1 



-Tr 



1 



N 



Ei-H E 2 -H c 6(#i - £ 2 ) 3/2 



const 



by/Ei 



(47) 



However, this expression evidently contradicts to the simple estimate of any of the diagrams of fig. 10. 

Thus, we have to choose between the two scenarios. First, as it is seen from the naive estimates the ~ 1/6 corrections 
to the correlation function may be much more singular than it is expected from the universality (|47| ) . In this case the 
~ 1/6 contributions blow up at AE much larger than the effective energy level spacing AE e ff ([13[) which will be the 
indication of some new physics at the intermediate energies AE e f f <C AE <C V . 

In the second scenario the universality (|^,|^]) do takes place and all the additional singularities cancel each other 
in the sum of different diagrams. 

In fact the aim of this section is to demonstrate explicitly that at least the ~ 1/6 corrections to the correlation 
function do not violate the universality (42|4^) and all the additional singularities vanish after a huge cancellation 
between the diagrams of fig. 10. 

The analogous cancellation between the high order corrections to the correlation functions has been previously 
observed for usual N x N random matrices by J.Verbaarschot et al. [0J|]. 

Consider first the most singular diagrams of fig. 10a. The calculation of the corresponding contribution to the 
correlation function of two logarithms (R3) has much in common with the calculation of the ~ 1/6 2 correction to the 
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Green function at the edge (|27||29j) . The main difference is that now two of the double- link chains are accompanied 
by the factor (V 2 Gq(Ei)Go(E 2 )) p = e~ Ap , while the third chain is associated with the oscillating contribution 
(V 2 Gl(Ei^)) k = exp(±2ik(j) — k\) . It is convenient to brought together into one expression the both diagrams of 
fig. 10a 



*M. = W* E E 



N 



Pl,P2>0 — oo<fc<oo 
^ e 2icf>k — \(p 1 +p 2 + \k\) 



e 2i<pk~\(p 1 +p 2 + \k\) 

(2Tr) d ( Pl p 2 \k\) d / 2 



1 

Pi 



1 

P2 



1 

W\ 



-d/2 



b 2 ^ (2ir) d ( Pl p 2 + \k\( Pl + p 2 )) d / 2 



(48) 



Naively summation over pi and p 2 here gives the singularity ~ A _2+rf (and the additional factor A -2 will appear after 
differentiation with respect to E\ and E 2 which should be done in order to get the correlation function of two Green 
functions). However, at least this leading singularity should disappear after summation over k. For example, for d = 
one has ~^2 k e 2l< M- A l fe l ^ \ (even more trivial example of the same kind is XX - l)* = 0)- For d ^ some cancellation 
also should take place at least for large p\. 2 which, on the other hand, are responsible for the "naive" singularity of 
(fi"8|). Thus again one may see that large k or long chains on the diagram (continuum limit) turn out to be important. 
In order to treat this cancellation explicitly it is convenient to divide the contribution of fig. 10a into two parts 



S ( LL )a 

A 



A + B , 
N x ^ 

b 2 (2TT) d ^ 

Pi,P2>0 — oo<A:<oo 



E 



,2i<j>k-X(p 1 +p 2 + \k\) 



( PlP2 + \k\( Pl + p 2 )) d / 2 



Vt! \V2 + |fc|) rf / 2 P d 2 /2 (Pi + \k\y/\ 



(49) 



B 



2N 



b 2 (2n) d 



-Xp 



.P>0 



d/2 



E E 

p'>0 -oo<fc<oo 



s 2i0fe-A(p' + |fc|) 
( P '+ \k\) d / 2 



Here in the part A the contribution with large both pi and p 2 will be suppressed due to summation over k, while the 
contribution with say pi <§C p 2 will be suppressed due to simple cancellation of the two terms in the figure brackets in 
( [49l ) . In the second part B summations over p-s are factorized and only summation over p' suffers from the cancellation 
due to oscillations. 

In the similar way one may write down the contribution of naively less singular diagrams of fig. 10b 



SILL) 



2N 



b b 2 {2it) d 
2N 



6 2 (27r)° 



E £ 



E 

.P>0 



-Xp 



d/2 



-Xp 



/c>0 q>0 



nd/2 



d/2 



E E 

-oo<fc<oo -|fc|<p'<0 



e 2i<pk-X\k\ 
(p'+\k\) d / 2 



(50) 



The situation is further simplified if one combines this contribution with the most singular part B (S3) of the fig. 10a 



diagrams. After simple change of variables p' 

2N 



5 LL) 



B 



6 2 (2tt) 



E £ 



P ■ 

'Xp 



\k\ one gets 



d/2 



(51) 



E _ 

- oc < fe < oo q>0 



E ^ m 4f2 *M-Hq \k\)9(q \k\)} 



Here 0{x) = for x < and 9(x) = 1 for x > 0. Now the sum over k may be found exactly 

S(LL) b 



B — 



2N 
b 2 (2ir) d 



-Xp 



lp>0 



d/2 



cos(2q<fr)e 



-Xq 



q>a 



g d / 2 2sin 2 



We see that this contribution as a function of A turns out to be as singular as the leading order result 
suppressed like 1/6 and therefore should be neglected. 



(52) 



)@ but is 



15 



Thus, let us consider the only surviving contribution A from (fi"9|). For simplicity consider the band matrices only 
(d = 1). After differentiation with respect to E\ and E% one gets 

5(Tr Gl TrG 2 ) = g2) £ JW-Mp^MM) {pi + p2 + \ k{ f (53) 

f 1 1 1 1 



[ \Jv\V2 + \k\(pi +p 2 ) \Jv\ip2 + |fc|) v^abi + l fc l) J 

As we will see both pi and j>2 in this sum effectively turns out to be large. Therefore, in order not to get exponentially 
small result one has to consider the singular in k contributions in the sum. These singularities naturally appear due 
to \k\ in (|53|). The following simple identity shows how one may utilize this ~ \k\ behavior: 



Y\k\e 2i * k ~ x W =Y\k\e 2i * k - x Mf(k) = 1- , (54) 

^— ' ^— ' 2 sin A 



where f(k) is any smooth and slow function of k and /(0) = 1. Taking into account that 2Vsin<f> = V AV 2 — E 2 one 
finds from j5^ , |54] ) 

d(jrG 1 J ? G 2 )- 62(4y2 _ ^ 2)2 2^ 27V _e — 6 2 4F 2 — i? 2 A£" 2 (55j 

Finally, the smoothed density-density correlation function for band matrices takes the form 



l p(Eh) p{E 2 ) 2 \\Ei_-Ez\J \ \\Et-E- 

where AE eff = N^irlpjE)]- 1 and I = b 2 {l-E 2 /AV 2 ) in accordance with @@||). For convenience we have added 
one factor (27r) _1 into the definition of AE e ff compared to (filf). 

To conclude this section let us remind again that we consider only the smoothed correlation functions. If one would 
like to compare our equation (^6|) with the result of numerical matrix diagonalization, the "experimental" result 
should be averaged with some smooth weight function. For example it may be 

/2 
K exp (AE + K x)e~ x -= , (57) 

where AE eff < k < AE. 



VI. SPATIALLY NONHOMOGENEOUS EXAMPLES 



The quantities which we have tried to calculate up to now - the density of states and density-density correlation 
function are generally considered for the usual N x N random matrices. In this section we would like to consider the 
two quantities which are specific for band matrices and never appear for the usual N x N ones. 
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n 2 

FIG. 11. Spatially nonhomogeneous correlation function p(Ei,i)p(E2,j). 

The first example will be the correlation function of local density of states p(E, i) for different energies and different 
vector indices (ld-lattice sites) i and j. The corresponding dressed Feynman diagram is shown in the fig. 11. Again 
one should calculate hrst the log-log correlation function. After differentiation with respect to E\ and E2 and taking 
the imaginary part of the Green functions the correlation function takes the form 



1 



p(E 1 ,i)p(E 2 ,j) c = -^G (E 1 )Go(E 2 ) 



Re 



E E 

fel,2>0 ni+712>0 



TT 2 b 2 

(Go (E 1 ) V) 2fcl (Go (E 2 )V) 2k - (Go (E 1 )G (E 2 )V 2 )^ +" 2 



(58) 



a/ (fci + fc 2 )(ni + n 2 ) + n\n 2 



1 

2~/7 



1 



exp 



2b 2 k 1 + k 2 



n 1 +n 2 



Here neither n\ nor n 2 could be negative. This equation is further simplified if one takes into account that effectively 
&i,2 "1,2- Therefore the k\ %2 may be neglected in the exponent and in the square root in the denominator and the 
sum over k\ and k 2 reduces to the simple geometrical progression. Finally, the summations over n\ and n 2 factorizc 
and the correlation function takes rather simple form 



p(E u i)p(E 2 ,j) c 



n 2 b 2 (AV 2 - E 2 ) 



Re 



E 



'27m 



exp 



{i-3? 
2nb 2 



(59) 



Here A is defined by (B8J). One may easily investigate for example small and large i — j limits of this expression. In 



terms of universal variables AE e f / and I (| 



p{E 1 ,i)p(E 2 ,j) c 
(pjE)/N) 2 



1 

Re 



correlation function 
2 



91) takes the form 



i-3 



x 



(60) 



exp 



1 Ei — E 2 



U AE, 



eff 



I 




Here Im(Ei — E 2 ) > and the integral should be squared before taking the real part. We have divided p(E) by N 
in the left hand side of (|6(]) in order to get rid of the physically trivial iV-dependance in the right hand side also. 
Finally, the integration over y in (|60| ) may be done explicitly, which leads to 



p(E 1 ,i)p(E 2 ,j) 
(p(E)/N) 2 



= -Re 
2 



iAE. 



Ei — E 2 



exp 




i — j\ sin(t)e 



(61) 
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where t = \i — j\/lyJ\E\ — E2\/ AE e f f . In particular one may easily examine that after summation over i and j the 
equation ([n]) leads to the usual correlation function (fiB]). 

Another interesting objects which may be considered easily within our technic are the finite size band matrices. 
Naturally the most interesting case is b <C A < b 2 . Below we describe analytically the crossover from the band matrix 
regime (^0|) to Wigner-Dyson regime in the asymptotics of smoothed density-density correlation function for finite 
size band matrices. 

Consider the periodic A x A band matrices. The statistical properties of this Gaussian ensemble are again defined 
by the second moment (JIJ), but now the function F takes the form 

+00 

F(*J)= E fii-3+nN) . (62) 



n— — 00 



Here f(k) = f(\k\) vanishes for k > b just like F in |IJ). The parameters V and b (strength of the interaction and 
width of the band) are now defined as 

N +00 2 f t \ 

7 2 =^F(i,i) = ^/(i) , b 2 = ±y^± , (63) 

j=l -00 



which is the natural generalization of (|2j). Also the analog of the equation ( |15[ ) for two link chain has the form 

N +00 

*„+l(i)=X>(M)*«0')= E /(i-J^nO*) > (64) 
3=1 3=— 00 

The solution of this equation for sufficiently large n (and for b <C AT) reads 

The leading order spherical Green function ( |To|) d ue to (^) (see also discussion before the equation ((!(])) is not 
changed. Therefore, the trivial modification of (|40[) gives 

= ^ - Re g /I e ' Ap exp { - ^ (x ) 2 } < (66) 

where A again is defined by (j38|). This expression may be further simplified in two limiting cases. If A is small, namely 
A -C b/y/X, one may replace summation over /c by integration. In this case the equation ([36]) reduces to the usual 
Wigner-Dyson correlation function p(E\)p(E2) c ~ 1 / [E\ — E2) 2 . If A 3> b/y/X only fc = survive in the second sum 
and d66| ) coincides with the pure band matrix result ( |45| ) . 

Also it may be convenient to use the "physical" variables: localization length I = b 2 (l — E 2 /4V 2 ) and effective 
interlevel spacing AE e ff = N[2ivlp(E)]~ 1 (note that AE e ff does not depend on A because p(E) ~ A), as it was 
done in (Bq) and (|60|). Now one has instead of (|66]) 



(p(E)) 2 V8^N J ^ \ At AE eff V I J 2x2 

/ / A£' (f/ V /2 f 1 _ r £\ >n(r)+siu(r) 



2 A V I -Ei - £2 1 / V dry ch(r) - cos(r) 



where r = f VPl = WfAEeff • In particular for t < 1 or t > 1 this equation corresponds to the usual Wigner- 
Dyson or band matrix (pq) results. 
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VII. CONCLUSIONS 



Random matrix models are usually expected to describe some universal and very general features of complicated 
quantum systems. Therefore on the one hand, one and the same very simple model may be associated with a variety 
of physical systems. On the other hand, this model will be generally able to explain only simplified versions of 
a real complicated problem, say properties of only very small metallic grains. In this paper we have considered 
banded random matrices which at least formally seem to be much closer to the real physical systems. For example, 
realistic Hamiltonians in the shell model for complicated atom pl| and atomic nuclei p2|] were shown to have a 
banded structure. Also, being a good example of quasi-ld quantum system random band matrices are expected to 
depict adequately properties of electrons in thick wires (see also |23| there the mapping of the Hamiltonian for 

disordered wire onto random banded block-diagonal matrix was done explicitly). 

Technically our work was stimulated in part by the successful application of matrix models for calculation of 
partition function of 2<i-quantum gravity. Just like in 2d-gravity we have found a critical behaviour at the edge of 
spectrum for band matrices (Section 3) and for lattice Hamiltonian with random hopping (Section 4). Unfortunately 
it may be not so easy to find a physical system whose global spectral properties will be described by the random band 
matrices with their almost semicircle density of states. Nevertheless, the interest in investigation of the edge behavior 
and the tails of spectral density have been demonstrated in recent papers (24||ll| . 

In almost all applications of random matrix theory one is interested in local characteristics of spectrum such as 
correlations of very close or even neighbouring energy levels. In Section 5 we have calculated the asymptotics of two- 
point correlation function of the density of states ( ]4C|j45| ) which is in agreement with the result of ref. fis} ] for energy 
level correlations in disordered metals. Moreover, together with the hypothesis of universality of spectral correlations 
|L9| our result (|45| ) allowed to estimate the localization length (46|) for random band matrices and this calculation 
seems to be much less complicated than known from literature p ] . 

On the other hand, we are able to calculate only the asymptotics (plus corrections) of the correlation function which 
in principle should not necessarily manifest universality. The universal behaviour of J45| ) shows that there are only two 
different energy scales in the model the global width of the energy zone AE = AV and the effective inter-level splitting 
AE e f f ( |43| , |56| ) . The hypothesis of universality finds further support in the calculation of the first correction to the 
two-point correlation function (|5^) . To the best of my knowledge it is the first calculation of subleading corrections 
for quasi-l<i systems. However, the calculation of the correction also shows a serious drawback of our perturbative 
approach. The accurate result ( |55| , [56| ) was found only after the huge two-step cancellation. One may speculate that 
this is the price for working very far from the region of convergence of the initial series in H/E. Nevertheless, these 
cancellations show that it will be extremely difficult to reach the region E\ — E% ~ AE e ff in our approach. 

Finally, in the last section we have found the asymptotics of the local density of states two-point correlation function 
( |60| ) as well as the usual two-point correlator for finite size quasi-ld system (^). These relatively simple analytical 
calculations demonstrate again the usefulness of diagrammatic approach for investigation of such nontrivial systems 
as random band matrices. 
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